Predicting a Gapless Spin-1 Neutral Collective Mode branch for Graphite 
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Using the standard tight binding model of 2d graphite with short range electron repulsion, we 
find a gapless spin-1, neutral collective mode branch below the particle-hole continuum with energy 
vanishing linearly with momenta at the Y and K points in the BZ. This spin-1 mode has a wide 
energy dispersion, to ~ 2 eV and is not Landau damped. The 'Dirac cone spectrum' of electrons at 
the chemical potential of graphite generates our collective mode; so we call this 'spin-1 zero sound' 
of the 'Dirac sea'. Epithermal neutron scattering experiments, where graphite single crystals are 
often used as analyzers (an opportunity for 'self-analysis'!), and spin polarized electron energy loss 
spectroscopy (SPEELS) can be used to confirm and study our collective mode. 
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Graphite is an important system in condensed mat- 
ter science and technology; in carbon research its role 
is fundamental. Its electrical and magnetic properties 
have been investigated for decades both experimentally 
and theoretically jjj . It is one of the simplest of quasi two 
dimensional zero gap semiconductors/semi metals. Inter- 
calated graphites offer many phases of condensed matter 
including superconductivity. Other important systems 
such as Bucky balls, carbon nano tubes and some form 
of amorphous carbon derive many of their novel prop- 
erties from their underlying 'graphite character'. Any 
newer understanding of graphite is likely to have a wider 
impact. 

The aim of the present letter is to predict a simple 
but important property of graphite that calls for re- 
examination of some of the low energy electrical and 
magnetic properties of graphite. We find that graphite 
possesses a new, unsuspected gapless branch of a spin- 
1 and charge neutral collective mode. This branch lies 
below the electron-hole continuum (figure 1); its energy 
vanishes linearly with momenta as Ku> s f=a hvpq(l — aq 2 ) 
about three points {T,K K') in the BZ (figure 1). 

Since graphite interpolates metals and insulators, our 
collective mode can be viewed both from metallic and 
insulating stand point. In paramagnetic metals 'zero 
sound' is a Fermi surface collective mode (3) . The 'charge' 
as well as the 'spin' of a Fermi sea can undergo indepen- 
dent oscillations. The charge oscillation becomes a high 
energy branch, the plasmon, because of the long range 
coulomb interaction; plasmons in graphite has been stud- 
ied in great detail in the past Q. The electron-electron 
interactions in normal metals do not usually manage to 
develop a low energy spin collective mode branch because 
of the nature of the particle-hole spectrum. However, the 
particle-hole spectrum of 2d graphite with a 'window' (fig- 
ure 2) provides an unique opportunity for a spin-1 collec- 
tive mode branch to emerge in the entire BZ. From this 
point of view our spin-1 collective mode is a 'spin-1 zero 



sound'(SZS) of a 2+1 dimensional 'massless Dirac sea', 
rather than a 'Fermi sea'. 

From an insulator point of view our collective mode is 
a spin-triplet exciton branch. Triplet excitons are well 
known in insulators, semiconductors and pit bonded pla- 
nar organic molecules; however, they usually have a finite 
energy gap, except when there are magnetic instabilities. 

Our spin-1 collective mode may be thought of as a 
manifestation of Pauling's || RVB state of graphite: the 
spin-1 quanta is a delocalized triplet bond in a sea of res- 
onating singlets. The gapiessness makes it a 'long range 
RVB' rather than Pauling's short range RVB. Later we 
will present an argument to suggest that at low ener- 
gies the neutral spin-1 excitation might undergo quantum 
number fractionization into two spin-i spinous. 

Existence of our gapless spin-1 collective mode branch 
should influence the spin part of the magnetic suscepti- 
bility, rather than the orbital part, which for graphite 
is diamagnetic, large and anisotropic. Study of spin sus- 
ceptibility by ESR, NMR and inelastic neutron scattering 
are good probes to detect the low energy part of our col- 
lective modes over a limited energy up to ~ 50 meV. 
A recent observation of 'large internal fields' in oriented 
pyrolitic graphite by Kopelevich and collaborators @ , in 
their ESR studies could be due to our low energy spin-1 
collective modes around the T point in the BZ. The low 
energy collective modes also contribute to specific heat 
and thermal conductivity over a wide temperature range. 
Our mode could be probed over a large energy range, by 
epithermal neutrons and spin polarized electron energy 
loss spectroscopy (SPEELS) @. In view of a wide energy 
scale associated with the collective modes, probes such as 
two magnon Raman scattering, ARPES, STM and spin 
valves Q should also be tried. 

Importance of electron-electron interaction in graphite 
P.pi)[ and related systems [|l3]-[l5| has been realized re- 
cently and it has lead to several interesting studies and 
predictions. 2d cuprates with Dirac cone spectrum has 
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been studied in the context of AFM order in the Mott in- 
sulating RVB-flux phase, for spin-1 goldstone modes |ll]] 
and d-wave superconducting phases, for spin-1 collective 
modes 12 1. 

Real graphite is a layered semimetal - stacked layers 
of honeycomb lattice of carbon atoms. We have one 
p z orbital as the relevant valence orbital and one elec- 
tron per carbon atom. The pir bond produces a filled 
valence band and an empty conduction band with van- 
ishing band gaps at two K points in the BZ. The coupling 
between graphite layers is van der Waals like. However, 
a small 'coherent' interlayer hopping has been invoked 
to explain the presence of small electron and hole tubes 
and pockets, (with 10 -4 carriers per carbon atom, i.e., a 
Fermi energy tp ~ 100 — 20QK ), responsible for the semi 
metallic character of graphite. 

We start with a 2-dimensional Hubbard model for 
graphite, which captures the physics of low energy spin 
dynamics. The Hamiltonian is: 

H = -t { c i,<? c j,o + M + u n ^ n H 

Here t ~ 2.5 eV is the nearest neighbor hopping ma- 
trix element. While the bare atomic U is of the order of 
8 eV, the effective renormalized U can be of the order of 
3 — 4 eV. We will keep U as a parameter to be fixed by 
experiment. 

The dispersion relation for the pn -bands is: 




^/?>k x a kyd kyd 
4 cos — - — cos + 4 cos 2 (1) 



with vanishing gaps at the two K points in the BZ (figure 
1). The particle-hole continuum of excitations is shown 
in figure 2. The 'Dirac cone single particle spectrum' at 
the r and K points makes the particle hole continuum 
very different from that of a free Fermi gas, or systems 
with extended Fermi surface. In contrast to figure 3, the 
particle-hole spectrum of a 2d Fermi liquid, our spectrum 
has a 'window'. The 'window' is characteristic of a Id 
particle-hole spectrum. In the Hubbard model two par- 
ticles with opposite spins at a given site repel with an 
energy U. This means an attraction for up spin particle 
and down spin hole; or an attraction in the spin triplet 
channel for a particle- hole pair. A spin triplet particle- 
hole pair could form a bound state, provided there is suf- 
ficient phase space for the attractive scattering. We find 
one spin-1 bound state for every center of mass momen- 
tum of the particle- hole pair. In particular an effective Id 
character of phase space also makes the collective mode 
energy vanish linearly with momenta around the three 
points: T and K's. 




FIG. 1. (a) Honeycomb lattice; ai, 2 = §(V3, ±1) (b) the 
Brillouin Zone; bi,2 = 2^(;4=,±1) (c) Dirac cone spectrum 
at a K point. 

The collective mode that we are after are obtained as 
the poles of the particle- hole response function in the spin 
triplet channel. We will focus on the zero temperature 
case. The magnetic response function within the RPA 
(particle-hole ladder summation) is given by: 

x°(q,w) 



(2) 



1 - u(q)x°(q,w) 
For Hubbard type on site repulsion , u(q) = U and the 



free particle susceptibility is : 
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Here /'s are the Fermi distribution functions. We have 
evaluated the RPA response function numerically and 
found the collective mode branch in the entire BZ, below 
the particle-hole continuum. However, it is instructive 
to linearize the electron and hole dispersion for low ener- 
gies, a Dirac cone approximation ||, and get an analyti- 
cal handle. We linearize the dispersion around K and K' 
and replace the BZ by two circles of radii k c (figure 1): 

e k = ±v F I k | fork < k c (4) 

where vp = -^i and N is the number of unit cells. In 
our linearization scheme, in equation 4 the summation is 
over the two circular patches (figure lb). 

For a finite range of q and cj, Im X°( uj j q) can be eval- 
uated exactly M: 



Im x°(q,w) = 
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with a square root divergence at the edge of the particle- 
hole continuum in (w,q) space. This expression has 
the same form as density of states of a particle in ID 
(with energy measured from vfo)- Note that in fact 
Im x°(<7,o;) = itp q (uj), where p q (u>) is free particle-hole 
pair DOS for a fixed center of mass momentum q. That 
is, the particle-hole pair has a phase space for scatter- 
ing which is effectively one dimensional. Thus we have a 
particle-hole bound state in the spin triplet channel for 
arbitrarily small U. However, we also have a prefactor 
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that scales the density of states. This together with 
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FIG. 2. (a)particle-hole continuum with a 'window' for 
graphite, (b) S(q,u>) for q > q c (~ if) 
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FIG. 3. (a)particle-hole continuum without a 'window' for 
a 2d fermi gas. (b) 5(q, cj) for q < 2&f. 



the square root divergence of the density of states at 
the bottom of the particle-hole continuum gives us a 
bound state for every q as q — > 0, with the binding en- 
ergy vanishing as aq 3 , as shown below. The square root 
divergence has the following phase space interpretation. 
The constant energy (Huj) contour of a particle-hole pair 
of a given total momentum q defines an ellipse in k-space: 
w = Ui?(|k + q|+|k|). In our convention, the points 
on the ellipse denote the momentum co-ordinates of the 
electron of the electron-hole pair. As the energy of the 
particle-hole pair approaches the bottom of the contin- 
uum, i.e., e p -h — > vfQ, the minor axes of the ellipses 
become smaller and smaller and the elliptic contours de- 
generate into parallel line segments of effective length 
~ g3. The asymptotic equi-spacing of these line seg- 
ments leads to an effective one-dimensionality and the 
resulting square root divergence. 

According to (1), the collective mode in magnetic 
channel is the solution of: 

l-[/ X (q,^) = 

or equivalently,/m x°(q, uo) — and Re x°(q, w) = jj. 
The asymptotic expression for Re \° (q> <*0 is found to be 
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where z = Using the above expression we obtain 

the following dispersion relation for the collective mode: 
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qv F - E B (q) (5) 



as uj — > q —> 0. Here Es{q) is the binding energy of the 
particle- hole pair of momentum q around the T point. 
The binding energy around the K points is roughly half 
of this. 

We mentioned earlier that our collective mode is a 
'magnetic zero sound'. While magnetic zero sound are 
difficult to get in normal metals, graphite manages to 
get it in the entire BZ because of the window in the 
particle- hole spectrum (figure 2). 



Having established the existence of a gapless spin- 1 col- 
lective mode branch within Hubbard model and the RPA 
approximation, we will discuss whether the semi-metallic 
screened interaction of 3d stacked layers will affect our 
result. As mentioned earlier, in tight binding situation 
like ours, the spin physics is mostly captured by the short 
range part of the repulsion among the electrons. We have 
numerically studied the response function for a more real- 
istic intra layer interaction namely the screened coulomb 
interaction (including interlayer scattering between lay- 
ers separated by distance d) given by B 



v(u),q) = 
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and find that the collective mode survives with small 
quantitative modifications. 

Let us discuss life time effects, that is beyond RPA. A 
remarkable feature of our collective modes is that it never 
enters the particle-hole continuum. It does not suffer 
from Landau damping (resonant decay into particle-hole 
pair excitations). To this extent our collective modes 
are sharp and protected; higher order processes will pro- 
duce the usual life time broadening, particularly at the 
high energy end. However, in real graphite there are 
tiny electron and hole pockets in the BZ with a very 
small Fermi energy ~ 10 to 20 meV. This leads to 'Lan- 
dau damping' of low energy collective modes around the 
r and K points, but only in a small momentum region 
Afc ~ 2fcp ~ -eh-, where kp is the mean Fermi momen- 
turn of the electron and hole pockets. That is, only a few 
percent of the collective mode branch in the entire BZ is 
Landau damped. 

A small interlayer hopping between neighboring layers 
t± ~ 0.2 eV (<< 2.5 eV,the in plane hopping matrix 
element), has been always invoked in the band theory 
approaches to understand various magneto oscillation ex- 
periments and also c-axis transport in graphite. However, 
a strong renormalization of t± is possible, as anomalously 
large anisotropic resistivity ratio ~ 10 4 have been re- 
ported in some early experiments on graphite single crys- 
tals; a many body renormalization is also partly implied 
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by the existence of our spin-1 collective mode at low en- 
ergies. As the emergence of the small electron and hole 
pockets (cylinders) are due to interlayer hopping, inter- 
layer hopping affect the spin-1 collective modes only in 
a small window of energy and ~ 0.1 eV. For the same 
reason the collective modes do not have much dispersion 
along the c-axis. 

Within our RPA analysis the collective mode frequency 
becomes negative at the T point for U > U c ~ It. Be- 
cause there are two atoms per unit cell, this could be 
either an antiferromagnetic or ferromagnetic instability. 
Other studies Jl^,[l6| have indicated an AFM instability 
for U >U C ~ 2t. 

Now we discuss the experimental observability of spin- 
1 collective mode branch. The collective mode has a wide 
energy dispersion from to ~ 2 eV. The low energy 
to 0.05 eV part of the collective modes determines the 
nature of the spin susceptibility (equation 3) of graphite 
and leaves its signatures in NMR and ESR results. For 
higher energies we have to use other probes. 

Inelastic neutron scattering can be used to study the 
line shapes and dispersion of our spin-1 collective modes. 
However, epithermal neutrons in the energy range 0.1 eV 
to ~ 1 eV, rather than the cold and thermal, 0.2 to 
50 meV , neutrons is needed in our case, due to the large 
energy dispersion. The dynamic structure factor 5(q, uj) 
as measured by inelastic neutron scattering is obtained 
by using our calculated RPA expression for our magnetic 
response function using the relation: 



(1 
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- Im x(q,w) 
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(6) 



At the present moment one need not concentrate on the 
energy resolution and it will be good to focus on proving 
the existence of the spin-1 collective mode by neutron 
scattering experiments. As the single phonon density of 
states of graphite vanish for energies > 0.2 eV, one need 
not perform spin polarized neutron scattering in order to 
avoid single phonon peaks. 

Another probe for studying the spin-1 collective mode 
is the spin polarized electron energy loss spectroscopy 
(SPEELS); exchange interaction of the probing electron 
with the 7r-electrons of graphite can excite the spin-1 
collective mode. As the electron current and spin de- 
polarization essentially measures the magnetic response 
function, our calculation of x(q, a;) (equation 2) can be 
profitably used to interpret the experimental results. 

The square root divergence of density of states at the 
bottom edge of the particle-hole continuum tells us that 
the low energy spin physics is effectively one- dimensional. 
To that extent, in a final theory ; we may expect our spin- 
1 excitation to be a triplet bound state of 'two neutral 
spin— g spinons' rather than ' e + e~ electron-hole pairs'. 
Further, as the energy of the spin-1 quantum approaches 
zero the binding energy also approaches zero and the 



electron-hole bound state wave function becomes ellipti- 
cal, with diverging size. We may then view the low energy 
spin-1 quanta as a 'critically (loosely) bound' two spinon 
state, very much like the quantum number fractionization 
of the des Cloizeaux-Pearson spin-1 excitation in the Id 
spin-i antiferromagnetic Heisenberg model. Our result 
also suggest a non-linear sigma model and novel 2 + 1 
dimensional bosonization scheme for graphite ]T^ | . 

We find Q that our spin-1 collective mode survives 
in carbon nano-tubes in a modified fashion. Preliminary 
study shows that three dimensional scmimctals Bi, HgTe 
and a-Sn do not have spin-1 collective modes at low en- 
ergies, because of quadratic dispersion at the zero gap. 
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